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Numerical renormalization-group results on entropy of the anisotropic two-channel Kondo 

^ , model with the band-width cutoff {D) in the presence of a magnetic field (h) are obtained 

to determine crossover temperature from the non- Fermi liquid to Fermi liquid fixed point, ft 

(y-^ I is found that the crossover temperature (Tx) is given by Tx = tTk ~ D{AJ/Jav)^e~^^'^'^"' 

CO ' when (/i/Tk)^ <C r- ^ 1, where Tk, Jav and AJ are the Kondo temperature, the average and 

O ' . 

' difference of the exchange coupling constants, respectively. This result indicates that non-Fermi 

■ liquid behavior can be seen even if A J > Tk. Robust similarities of the crossover behavior in 

, the region around the non-Fermi liquid critical point to that of the two-impurity Kondo model 

^ , are also discussed. 

^ : 
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, §1. Introduction 

Non-Fermi-liquid (NFL) behavior observed in some heavy fermion compounds and High 
cuprates has stimulated intense studies for a variety of impurity models which have NFL quantum- 
critical points. ^"^^ While each of the models has its own energy scale Tk below which the ther- 
modynamic and transport properties are geverned by the NFL fixed point, relevant perturbation 
which is normaly present in real systems introduces another energy scale where crossover to the 
Fermi-liquid (FL) fixed point occurs. In other words, the NFL behavior can be observed only when 
Tx/Tk <^ 1. Since real systems would not be exactly on the critical points, it is important to know 
about the above condition that NFL behavior occurs. In this paper, we investigate the crossover 
energy scale Tx so that the condition Tx/Tk ^ 1 is specified in a microscopic sense by using the 

* E-mail: yotsu@eagle.mp.es.osaka-u.ac.jp 
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numerical renormalization-group (NRG) method^"^) for two types of the quantum impurity model, 
the two-channel Kondo model in the presence of a magnetic field and the two-impurity Kondo 
model, with a special emphasis on the effect of channel anisotropy. 

The non- Fermi- liquid fixed point of the multichannel Kondo model in the overscreened case, i.e., 
when the number of channels n is larger than twice the size of impurity spin S, is characterized 
by a degeneracy of the ground state which cannot be completely lifted by coupling to conduction 
electrons. ^'^-^ The residual entropy is given by ln[sin[(25 + l)^/(n + 2)]/sin[^/(n + 2)]].!°'") In 
temperature T below the Kondo temperature Tk, the specific heat coefficient 7 and magnetic 
susceptivility x are ~ 2^4/(n+2)-i^ ^.j^^ electrical resistivity p is const, it j'2/("+2) _io-i3) ^ magnetic 
field or channel anisotropy is relevant so that it brings about crossover to Fermi-liquid fixed points 
but 7, X and p keep the above scaling laws when <^ T <^ Tk. Fig.l presents a schematic 
phase diagram in the case of n = 2, S = 1/2 in which 7, x ~ and p ~ const, it -v/T 

in the region <ti T <^ Tk. The two-channel Kondo model has been regarded as one of the 
most probable candidates responsible for the NFL behavior observed in U alloys.'^) Recently, the 
multichannel Kondo model with channel anisotropy has been discussed in relation to singular effects 
of impurities in the region around the ferromagnetic quantum-critical point. ^'^^ 
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Fig. 1. A schematic phase diagram near a quantum-critical point (QCP) of the n = 2, 5 = 1/2 channel anisotropic 
Kondo modeL In the region Tx ^ T <C Tk, NFL behavior occurs due to the overscrening Kondo effect. 



The crossover behavior to the FL fixed point in the n = 2, S = 1/2 Kondo model with channel 
anisotropy and/or in the presence of a magnetic field has been investigated by a variety of meth- 

ods.15-22) It 

can be easily seen in the impurity entropy which decreases from In y/2 to around Tx 
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as temperature is lowered. The crossover temperature is proportional to the square of relevant 
perturbation as follows: 

Tx oc /l^ (1.1a) 
in the presense of a magnetic field H where h = hbH; 

oc (Ji - J2)^ (1.1b) 

in the channel anisotropic case where Ji and J2 are the exchange couplings for two channels, 
respectively.^^' On the other hand, the impurity free energy .Fimp has the following scaling form: 

Fnap(.r;T,h) = -Tf{r,T/TK,h/T), (1.2) 

where r is a dimcnsionless function of D, Ji and J2 which depends on a cutoff scheme and in general 
unknown. Since r represents channel anisotropy, r = if Ji = J2. From eqs. (1.1a) and (1.2), we 
can easily see that ~ h? /Ty, so that crossover from NFL to FL regime occurs if {h/Tx)'^ ^ 1 in 
the isotropic channel case.^' ^ ''^ In the case of channel anisotropy, however, eqs. (1.1b) and (1.2) 
tell us nothing about the condition that the NFL behavior can be seen. Our main purpose is to 
identify this condition analyzing temperature dependence of the entropy by the NRG method and 
it will be emphasized that crossover from the NFL to FL regime can occur even if | Ji — J2I > Tk. 

The two-impurity Kondo model is also known to have a NFL critical point, which arises from 
competition between the intcrimpurity antifcrroniagnctic interaction and single-channel Kondo 
effects. '^^ The nature of criticality of this system has been well understood.'^' The residual 
entropy is ln\/2 at the critical point just like the n = 2, S = 1/2 single impurity Kondo model 
which was discussed in the above. The staggered impurity susceptibility is | In T\ when 
Tx <C T <C Tk- Recently, it has been shown that a similar mechanism works in the model with the 
competition between the Kondo singlet and /^-crystalline-clectric-ficld singlet explaining the NFL 
behavior observed in Ri_^U^Ru2Si2 (R=Th, Y and La, x < O.O?).^'^'^^) 

The crossover temperature of the two-impurity Kondo model is known to be proportinal to 
{Kc — K)'^ where K is the intcrimpurity antiferromagnetic (RKKY) coupling constant and K^. is 
the critical coupling. We will also elucidate an explicit expression of the crossover temperature of 
this model and discuss similarities of the crossover behavior to that of the two-channel anisotropic 
Kondo model in a magnetic field. 

The organization of this paper is as follows. In §2, we discuss effects of channel anisotropy and 
a magnetic field in the two-channel Kondo model. In §3, we study the crossover behavior of the 
two-impurity Kondo model around the criticality. Finally, we summarize our results in §4. 
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§2. Anisotropic Two-Channel Kondo Model in a Magnetic Field 

2. 1 Model and numerical procedure 

In this section, we discuss crossover from the NFL to FL regime in the n = 2, 5 = 1/2 Kondo 
model with channel anisotropy and in a magnetic field described by the following Hamiltonian: 

H = Hk + Hi + Hz, (2.1) 

Hk= Y.^kcl^aCkma: (2-2) 

m=l,2 k,rj 

Hl= Y ^ Jmcl^^^aa'Ck'ma' ■ S, (2.3) 
m=l,2 k,k',a,a' 

Hz = -2hS^. (2.4) 

Here m (= 1, 2) is the index of two channels, a is the spin index, cj,^^ is the creation operator of the 
conduction electron, S is the impurity spin and a is the vector formed by the Pauli matrices. Hz 
represents the Zeeman term where the g'-factor is taken as 2. When Ji=J2 and h=0, one obtains 
the isotropic two-channel Kondo Hamiltonian leading to the NFL stable fixed point. The case of 
J1/J2 corresponds to the anisotropic two-channel Kondo Hamiltonian whose stable fixed point is 
that of FL. 

To analyze the Hamiltonian eq. (2.1) by the NRG method,^"^) we transform eqs. (2.2) and (2.3) 
into 

m=l,2n=0 T 
ft 



^{frn,nafi'hn+la + f m,n+la f'rn,na-) 1 (2-5) 



Hi = ^ ^ ^— Y'^"^fL,0a^<7<7'fm,0a' ' S, (2.6) 

m=l,2 a,a' 

where fm,n {fm n) annihilation (creation) operator of the conduction electron in the Wannicr 
representation whose extent /c^^A"/^. A is a logarithmic discretization parameter in NRG calcu- 
lation, D is half the bandwidth of the conduction electrons. Following an usual procedure of the 
NRG method, we define iJjv as 

Af-l 

H^ = A(^-i)/^ [E E E A-"/^(4„./.,„+i. + h.c) 

m a n=0 

+ fli,Oa^T,a'fm,Oa' ' S . (2.7) 

Then, the Hamiltonian (2.7) satisfies the recursion relation 

Hn+I = A^/^Hn + Y^J'ii,Nafm,N +lc7 + fm^N+lafm,Na)- (2.8) 

m,(T 

By repeated use of the recursive procedure, we can compute the eigen states and energy levels of 
Hamiltonians (Hn)- 
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Once the energy levels of iJjv are obtained, the entropy can be calculated by the relation 



S = PnHHn) - Fn), (2.9) 
in each AT-step. Here {H^) and are defined as 

Fn = --^ In^AT = --^ ln('IVe-''^-f^^), (2.10) 
Pjv Pat 

where /?jv = A^^^^^^/^/T. By setting T ~ Tjv = A~(-^~^)/^, we can obtain temperature dependence 
of the entropy with a good accuracy in NRG calculation. 

In our calculation, we set /3^^ = 0.5 and A = 3 and about 600 total states are retained at each 
iteration. Hereafter, we take the unit of energy as (1 + A'~^)D/2. 

2.2 Crossover temperature 

In this subsection, we elucidate an explicit expression of the crossover temperature of the 
model, eq. (2.1), from NRG results on the impurity entropy SirapiT). As mentioned in §1, if 
0<Tx/TK<d, SiinpiT) of this model decreases from 1/2 • In 2 to around as T is lowered. So 
can be defined as temperature at which the entropy takes a value of S'inip(?x)=l/4-ln2, the middle 
point between NFL and FL fixed points. (For examples, see Fig. 2 below.) 

Firstly, we investigate an effect of channel anisotropy on Tx for h = 0. In this case, we can assume 
that Tx is described by a function of two parameters, i.e., Jav = (Ji + ■h)/'^ and A J = Ji — J2. 
Our procedure to determine Tx is illustrated in Fig. 2 (a) and (b) which present NRG results on 
Simp{T) for Jav = 0.1, AJ = 1.0 X 10-4 and for Jav = 10, AJ = 1.0 x 10"^ respectively It is 
noted that Tx/Tk <C 1 while AJ^Tk in both (a) and (b). These examples lead us to the fact 
that Tx is not scaled as Tx ~ A J^/Tk. (If Tx were ~ A J^/Tk, AJ should be much less than Tk 
for Tx/Tk ^1.) Actually, we have tested a bunch of parameter sets and found a function which 
reproduces the behavior of Tx as follows: 

Tx = a X (AJ/Jav)'e-V'^-, (2.11) 

where a is a numerical constant which is taken to be a = 2.95 within the accuracy of our calculation, 
and 

AJ = ^aAJ, (2.12) 



Ja,v — ^A«^av) (2-13) 

1 1 + A- ^ 
21- A- 



^A= o^^^r^ilnA. (2.14) 
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Fig. 2. NRG results on the entropy S'imp(T) of the anisotopic two-channel Kondo model for (a) Jav =0.1, A J = 
1.0 X 10~^ and for (b) Jav = 10, AJ = 1.0 x 10~^. In each case, crossover from NFL to FL regime can be clearly 
seen, i.e. rx/rK<Cl, and then Tx is obtained by the relation 5'imp(rx) — 1/4 • In 2. 



Here A\ is a revision factor of the density of states in NRG calculation. Fig. 3 shows that is well 
fitted by eq. (2.11). From Fig. 3 (a) in which Jav=0.3 is fixed, Tx can be seen to be proportional to 
A J^, which is consistent with the results of Pang and Cox.^^) On the other hand, from Fig. 3 (b) 
in which A J=1.0 x 10"'^ is fixed, Jav-dependence of Tx can be seen to be equivalent to that of the 
right hand side of eq. (2.11). In the weak-coupling case of Jav^O.4, Tx is an increasing function of 
Jav as can be seen in Fig. 3 (b) while Tk is well known to be an increasing function ~ Jsm^~^^'^'^^ 
in the case of n = 2 by the two-loop perturbative renormalization group theory. It is noted that if 
Tx were ~ AJ^/Tk, Tx should be a decreasing function of Jav in the weak-coupling regime. (See 
also Fig. 4 (b) for comparison.) 

We have checked that eq. (2.11) can be adopted in a wide parameter region of Jav for the n = 2, 
5=1/2 Kondo model with channel anisotropy. In the weak-coupling regime, from eq. (2.11) and 
Tk ~ JavC"^/'^'"^, the condition that crossover from NFL to FL regime occurs, i.e. Tx/Tk^1, is 
given explicitly by AJ^/Jay<Cl. This result indicates that the NFL behavior in the two-channel 
Kondo model is fairly robust against channel anisotropy. 

Next, we investigate an effect of a magnetic field in the isotropic channel case of Ji= J2= J. In this 
case, as mentioned in §1, Tx is ~/i^/Tk.^'"^^'^'^^ As shown in Fig. 4 (a) and (b), it can be actually 
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Fig. 3. Crossover temperature (Tx) of the two-channel Kondo model with channel anisotropy (AJ). (a) Tx versus 
AJ with Jav=0.3 fixed, (b) Tx versus Jav with AJ=1.0 x 10"^ fixed. Tx is well fitted by eq.(2.11). 

checked from our NRG results on the entropy that can be well fitted by a function of 

= /3 X /iVlk, (2.15) 

where /3 is a numerical constant of /? = 4.45. Here Tk is defined in terms of the dynamical impurity 
susceptibility x{^) by 

Tk oc lim lim [Imyfo;)!^^ , (2.16) 

for the isotropic channel model in the absence of a magnetic field^^"^^^ in which the value of Tk is 
normalized to coincide with Je^^l"^ (=^k) ^=0.1. Eq. (2.15) is valid in a wide region from weak 
to strong coupling. This suggests that eq. (2.16) is a good definition of Tk- In Fig. 5, we show 
J-dependence of Tk, which is defined by eq. (2.16), and of T^, which is the two-loop perturbative 
renormalization group result. While T^ is a monotonously increasing function of J, Tk exhibits a 
broad peak at or around J~1.3. The discrepancy between Tk and T^ indicates that the two-loop 
perturbative renormalization group result on the Kondo temperature is not sufficient for J^O.5. 

2. 3 Scaling function for entropy in the crossover regime 

We have obtained an explicit expression of Tx, eq. (2.11), in the conventional band-width cutoff 
scheme for the n = 2, = 1/2 Kondo model with channel anisotropy. On the other hand, the 
thermodynamic Bethe ansatz (TBA) equations for the anisotropic multichannel Kondo model was 
obtained by Andrei and Jerez. While their cutoff scheme differs substantially from ours, and so 



7 



(X 10") 
300 



. (b) 



h = 10" 
■ Tx (NRG results) 
4.45x/i'/rK 




Fig. 4. Crossover temperature (Tx) of the two-channel Kondo model in the presense of a magnetic field {h). (a) Tx 
versus h with J—Q.Z fixed, (b) Tx versus J with /i=1.0 x 10~* fixed. Tx is well fitted by eq. (2.15). 
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Fig. 5. J-dependence of Tk defined by eq. (2.16) compared with T^=Je . 



does their expression of (7k)) their n = 2, 5" = 1/2 model must be in the same universahty class 
as ours. It is thus possible to make a comparison of our NRG results with the TBA result by use of 
our expression of (if necessary, Tk). In this subsection, we concentrate on the crossover regime 
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in which O^S'iinp(T')^ ln\/2, T^Tk, and exphcitly show that our NRG results on the entropy can 
be fitted by a scaUng function oiT/T^. 

First, we derive the scahng function in the crossover regime analytically from the TBA equations. 
In the scaling limit, the impurity part of the free energy is given by^^^ 



^ imp 



^^>- 27ry_ooCosh[e + ln(r/rK)]''^' ^ 



where r/i(^) is determined by the TBA equations; 

lnr/i(0 = -2re« + G* ln[l+ 7/2(0], (2-18) 
lnr?2(0 = -e« + G * ln[l + 771(6] + G * ln[l + r?3(0], (2-19) 
lnry,(0 = G*ln[l + r?,_i(e)]+G*ln[l + ??,+i(0], j > 3, (2.20) 

with a boundary condition 

Inr7,(£) 2h 

lim — = — . (2.21) 

j— S-OO J T 

Here G*(^(0 = /^oo G{i- and ^(^-^0 = l/27rcosh(e-e')- For T « Tk, the impurity 
free energy, eq. (2.17), is determined by the ^— >+oo asymptotic form of rji. In the limit of ^— >+oo, 
r]j is given by 

..(+oo) = ^^^^%lMl-l, j>2. (2.22) 

Firstly, we consider the case of T<^h<^TK- By eq. (2.22), 773(+oo) ~ exp(2^/r) and its correc- 
tions due to ^ dependence can be neglected, so that eq. (2.19) leads to 

\nr,2{0 = -e^ + h/T + f{r,e-^,h/T), (2-23) 

where / is defined by eq. (1.2). In eq. (2.23), the last term can be seen to compare negligibly with 
h/T if we note that this term is ~Fiinp(r)/r in a region around ^~ln(rK/r). Then, after some 
manipulations of eqs. (2.17), (2.18) and (2.23), we obtain 

F (T\ ^ ln[l + exp(-(Vr)e(C))] 
27r J -00 cosh[(; + ln(/i/rK)J 

where e{C,) is given by 




(2.25) 



~ e"^ (2r - e-2^ /27r) , for e""^ ~ H/Tk < 1. (2.26) 

The integral in eq. (2.24) has dominant contributions from a region around C ~ ^^{T^/h). When 
r S> (/i/Tk)^, the second term can be neglected in eq. (2.26) so that the impurity free energy can 
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be written in the following scaling form: 

i^imp(r) = -Tf{T/T^), (2.27) 
where the crossove temperature is given by = vTk and / is a scaling function in the crossover 



regime; 



f( N r In l + exp -2exp 0^ ^^ 

= loo 2.cosh[g + In(x)] ('-^'^ 

= — In^^ ^--Inrf— + +^ln27r, (2.29) 

irX TTX TTX V TTX 2 / 2 

where F is the gamma function. Then, the impurity entropy, S'imp=— ^^imp/^T", is calculated as 

= W/T7}''' (wy 2) " wJt3 

where -0 is the digamma function. On the other hand, when r <C {h/TyC)'^, the first term can 
be neglected in eq.(2.26). In this case, if the crossover temperature is regarded as Tx=h'^ /4ttTk, 
T-dcpcndent contributions in the impurity free energy is given by —Tf(T/Ty^) and therefore eq. 
(2.30) is also available for the entropy. 

Secondly, we consider the case of /7.^T<CTk. By eqs. (2.19) and (2.22), \nri2{Cj——c^^ so that 
lnr/i(^)~— 2rc^. Hence eq. (2.30) is available by Tx=rTx- It is noted that if r <C (/i/Tk)^, 
Tx=rTK<^h^T so that the entropy is ~ln ^/2, which is not sensitive to the value of Tx- 

In summary, a scaling function for entropy in the crossover regime is given by eq. (2.30) where 
Tx=rTK for r>(V7k)^ and Tx=/iV4vrTK for r<(/i/7k)^. For eq. (2.30), 5'imp(Tx)=0.99x 1/4-ln 2 
so that Tx is about the same as we defined in the former subsection. (We must apologize for the 
normalization schemes on Tk which are different between the former and this subsections.) 

Second, we compare our NRG results on entropy with eq. (2.30). Figures 6 and 7 present the 
results on entropy versus temperature scaled by T^, expressions of which are eq. (2.11) and eq. 
(2.15), in the anisotropic channel case and in the presence of a magnetic field, respectively. In both 
cases, the T dependence of entropy can be fitted quite well in the crossover regime by the scaling 
function, eq. (2.30), derived from the TBA equations. 



§3. Two-Impurity Kondo Model 

3.1 Model 

In this section, we study the two-impurity Kondo model whose Hamiltonian is given by two 
independent Kondo models coupled by interimpurity exchange interaction (K) as follows: 

H = HK + Hi + Hi^, (3.1) 



10 




10"' 10"' 10"' 10° 10' 10' 10' lO* 10' 
TIT^ 



Fig. 6. Temperature dependences of entropy obtained by the NRG method for the anisotropic two-channel Kondo 
model. In the crossover regime, by use of an expression of Tx, eq. (2.11), they can be fitted quite well by a scaling 
function, eq. (2.30) derived from the TBA equations. 




Fig. 7. Temperature dependence of entropy obtained by the NRG method for the two-cannel Kondo model in the 
presence of a magnetic field. In the crossover regime, by use of an expression of Tx, eq. (2.15), they can be fitted 
quite well by a scaling function, eq. (2.30) derived from the TBA equations. 
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Hk= Jl^kci^a^^kma, (3-2) 

m=l,2 k,(T 

Hi = Jm ^ cl.^^aaa'Ck'ma' ' Sm, (3-3) 

m=l,2 k,k',(T,a' 

= KSi ■ S2. (3.4) 

The notations are the same as before. In the original form of the two-impurity Kondo model studied 
by Jones et al.,^^ the Kondo exchange coupling between the conduction electrons and localized spin 
with both the same and different symmetry (specified by m) are taken into account, because they 
considered the case of the two impurities on different sites. In the present study, we neglect the 
Kondo exchange coupling between different symmetries, for simplicity. It is however noted that our 
model, eq. (3.1), itself can be considered to be a local-moment version of the two-orbital Anderson 
model with crystalline-electric-field (CEF) effect in a certain pseudospin representation.^ '''^^^ 

The model, eq. (3.1) has a NFL fixed point at K=Kc{Ji, J2) for a set of Ji and J2 due to com- 
petition between the single channel Kondo effects and the interimpurity antiferromagnetic interac- 
tion.^' In this section, we concentrate on the region around critical points of K=Kc{J, J). 
The definition of is the same as that of the former section; Simp (^x)= 1/4 • ln2. 

3.2 Crossover temperature in the region around critical points of K=Kc{J,J) 

Firstly, we consider cases in which K is displaced slightly away from a critical value of 
Kc=Kc{J,J) with Ji=J2=J fixed. It has already been discussed by using the conformal field 
theory that Tx is proportional to AK'^/Tk,'^^^ where AK=Kc — K. Here we define Tk by 

Tk = lim lim 1/Imxs(^<^), (3.5) 

ijj — ^0 T — ^0 

where Xs is the dynamical impurity susceptibility for — ^2. We have examined various parameter 
sets and confirmed that the crossover temperature Tx is given as 

Tx = 7 X Ai^V^k, (3.6) 

where 7 is a numerical constant of 7=0.50 within the accuracy of our calculation. In fact, as shown 
in Fig. 8, the temperature dependence of entropy can be scaled by Tx defined by eq.(3.6) quite 
well. 

Next, we consider cases in which J2 is displaced slightly away from J2=J with Ji = J and 
K=Kc=Kc{J,J) fixed. In Fig. 9, we show the behavior of a crossover for Ji=0.3, J2=0.297, 
and i^=i^c=0. 195416, compared with the case of the two-channel Kondo model with Ji=0.3 and 
J2=0.297. The isotropic channel cases in both models are also shown in this figure. In both cases, 
the small perturbation (Ji — J2) does not affect the high temperature behavior. Namely, it does 
not change Tk- The temperature dependence of entropy in the anisotropic channel case is quite 
similar in both cases in the crossover regime. Hence Tx for the two-impurity model is also given by 
eq. (2.11). 
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Fig. 8. Temperature dependence of entropy for the two-impurity Kondo model when K is displaced slightly away 
from Kc. Jfc^O. 564723 for Ji = J2=0.4. 7<'c=2. 655538 for Ji = J2=0.7. The temperature is scaled by T^, which is 
determined by eqs. (3.5) and (3.6). 
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Fig. 9. Crossover behavior for the two-impurity Kondo model (Ji=0.3, 72=0.297, and _ftr=ii'c=0. 195416), compared 
with the case of the anisotropic two-channel Kondo model with Ji=0.3 and J2=0.297. 
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§4. Summary 

In this paper, we have investigated effects of relevant perturbations, which bring about crossover 
from the NFL to FL regime, in the regions around quantum critical points of the two-channel and 
two-impurity Kondo models. 

For the two-channel Kondo model in which Tk is defined by Tk oc limj^^o limr-»o[Imx(^)]~^ (in 
the weak-coupling regime, TK—Ja,ye~^^'^'^"'), we have derived numerically exact expressions of the 
crossover temperature from our NRG results on the entropy as follows: (1) in the anisotropic 
channel case, if Tx/Tk^I (in the weak-coupling regime, AJ^/jfy<Cl ), crossover from the NFL 
to FL regime occurs around Tx which is given by Tx=ax{AJ / Ja,y)'^e~'^^'^'^'^ with a=2.95, where 
Jav and A J are the average and difference of the exchange couplings Ji and J2, respectively (a 
swung dash means A\ times, Aa being a revision factor of the order of 1); (2) in the presence of 
a magnetic field, if {h/Tx)^'^!, crossover from the NFL to FL regime occurs around which is 
given by T^=P x /tV^K, where /3=4.45 in our normalization procedure on Tk- 

For the two-impurity Kondo model which has a NFL critical point at K=Kc{Ji,J2) for a set 
of the Kondo exchange couplings Ji and J2, we have obtained NRG results on the entropy in the 
region around critical points of K=Kc{J, J) as follows: (1) when J2 is displaced slightly away from 
J2=J with Ji= J and K=Kc{J, J) fixed, crossover to a FL fixed point occurs around which is 
given by that of the two-channel Kondo model with the same set of Ji and J2 in the absence of 
a magnetic field; (2) when interimpurity (RKKY) exchange coupling K is displaced slightly away 
from the critical value of Kc{J, J) with Ji=J2=J fixed, crossover to a FL fixed point occurs around 
Tx which is given by Tx=^ x {Kf. — K)"^ /Ty, with 7=0.50 where rK=lim^^-->olimr->o I™Xs('^)- 

In the present work, all the results on Tx can be adopted in a wide region of J (Jav) and all the 
temperature dependences of entropy in the crossover regime can be fitted quite well by a scaling 
function, eq. (2.30), derived from the TBA equations of the two-channel Kondo model with channel 
anisotropy and in the presence of a magnetic field. It is emphasized that the NFL behavior in the 
two channel Kondo model is fairly robust against channel anisotropy so that crossover from the 
NFL to FL regime can be seen even if AJ>rK. 

Acknowledgments 

We wish to thank C. M. Varma for valuable advice. We also would like to thank I. Affleck, 
D. L. Cox, H. Kohno, H. Kusunose and K. Miyake for helpful comments. H. M. is supported by 
Reserch Fellowship of Japan Society for the Promotion of Science for Young Scientists. This work is 
supported in part by the Grant-in- Aid for COE Research (10CE2004) of the Ministry of Education, 
Science, Sports and Culture. 

1) P. Nozicrcs and A. Blandin: J. Phys. (Pairs) 41 (1980) 193. 

2) D. L. Cox: Phys. Rev. Lett. 59 (1987) 1240. 



14 



B. A. Jones, C. M. Varma and J. W. Wilkins: Phys. Rev. Lett. 61 (1988) 125. 

I.E. Porakis, CM. Varma and A.E. Ruckenstein: Phys. Rev. Lett. 70 (1993) 3467 

T. Giamarchi, CM. Varma, A.E. Ruckenstein and P. Nozieres: Phys. Rev. Lett. 70 (1993) 3967. 

K. G. Wilson: Rev. Mod. Phys. 47 (1975) 773. 

H. R. Krishna-Murthy, J. W. Wilkins and K. G. Wilson: Phys. Rev. B 21 (1980) 1003; 21 (1980) 1044. 
O. Sakai, Y. Shimizu and T. Kasuya: J. Phys. Soc. Jpn. 58 (1989) 3666. 

For example, see D. L. Cox and A. Zawadowski: Exotic Kondo Effects in Metals: Magnetic ions in a crystalline 

electric field and tunnelling centres (Taylor & Francis, London, 1999). 
N. Andrei and A. Dcstri: Phys. Rev. Lett. 52 (1984) 364. 

A. M.Tsvelick, J. Phys C18 (1985) 159. 

L Affleck and A.W.W. Ludwig: Nucl. Phys. B360 (1991) 641. 
L Affleck and A.W.W. Ludwig, Phys. Rev. B 48,7297 (1993). 
H. Maebashi, K. Miyake and C M. Varma: preprint, cond-mat/0109276. 

P. D. Sacramento and P. Schlottmarm: Phys. Rev. B 43(1991) 13294. 
H. B. Pang and D. L. Cox: Phys. Rev. B 44(1991) 9454. 

L Affleck, A. W.W. Ludwig, H. B. Pang and D. L. Cox: Phys. Rev. B 45 (1992) 7918. 

M. Fabrizio and A. O. Gogolin: Phys. Rev. B 50 (1994) 17732. 

M. Fabrizio, A. O. Gogolin and P. Nozieres: Phys. Rev. Lett. 74 (1995) 4503. 

M. Fabrizio, A. O. Gogolin and P. Nozieres: Phys. Rev. B 51 (1995) 16088. 

N. Andrei and A. Jerez: Phys. Rev. Lett. 74 (1995) 4507. 

P. Coleman and A. J. Schofield: Phys. Rev. Lett. 75(1995) 2184 . 

B. A. Jones and C. M. Varma: Phys. Rev. B 40 324 (1989). 
O. Sakai and Y. Shimizu: J. Phys. Soc. Jpn. 61 (1992) 2333. 

L Affleck and A. W. W. Ludwig: Phys. Rev. Lett. 68 (1992) 1046. 

L Affleck, A. W. W. Ludwig and B. A. Jones: Phys. Rev. B 52 (1995) 9528. 

S. Yotsuhashi, H. Kusunose and K. Miyake: J. Phys. Soc. Jpn. 70 (2001) 186. 

S. Yotsuhashi, K. Miyake and H. Kusunose: J. Phys. Soc. Jpn. 71 (2002) 389. 

D. L. Cox: Physica C 153-155 (1988) 1642. 

V. J. Emery and S. Kivelson: Phys. Rev. B 46 (1992) 10812. 

V. J. Emery and S. Kivelson: Phys. Rev. Lett. 71 (1993) 3701. 

Y. Kuramoto: Transport and Thermal Properties of f-Electron Systems, eds. G. Oomi, H. Fujii and T. Fujita, 
(Plenum Press, NewYork, 1993) p. 237; Y. Shimizu, O. Sakai, Y. Kuramoto: Physica B 206&;207 (1995) 135. 



15 



